Virial Masses of Black Holes from Single Epoch Spectra of AGN 



Brandon C. Kelly and Jill Bechtold 
bkellySas . arizona . edu , jbechtoldOas . arizona . edu 

Steward Observatory, University of Arizona, 933 N Cherry Ave, Tucson, AZ 85721 

ABSTRACT 

We describe the general problem of estimating black hole masses of AGN by calculating the 
conditional probability distribution of Mbh given some set of observables. Special attention is 
given to the case where one uses the AGN continuum luminosity and emission line widths to 
estimate Mbh, and we outline how to set up the conditional probability distribution of Mbh 
given the observed luminosity, line width, and redshift. We show how to combine the broad 
line estimates of Mbh with information from an intrinsic correlation between Mbh and L, 
and from the intrinsic distribution of Mbh, in a manner that improves the estimates of Mbh- 
Simulation was used to assess how the distribution of Mbh inferred from the broad line mass 
estimates differs from the intrinsic distribution, and we find that this can lead to an inferred 
distribution that is too broad. We use these results and a sample of 25 sources that have recent 
reverberation mapping estimates of AGN black hole masses to investigate the effectiveness of 
using the C IV emission line to estimate Mbh and to indirectly probe the C IV region size- 
luminosity {R-L) relationship. A linear regression of logLA(1549A) on logMnn found that 
1/1549 oc M^-^^^°-2^ A linear regression also found that Mbh oc L°^l^° °'' FWHM^jy, implying 
a C IV R-L relationship of the form Rciv oc L^glg^" "^. Including the C IV hue FWHM 
resulted in a reduction of a factor of ~ 1/3 in the error in the estimates of Mbh over simply 
using the continuum luminosity, statistically justifying its use. We estimated Mbh from both C 
IV and II/3 for a sample of 100 sources, including new spectra of 29 quasars. We find that the 

1/2 

two emission lines give consistent estimates if one assumes R oc Ljjy for both lines. 

Subject headings: galaxies: active — line: profiles — methods: data analysis — methods: statistical 
— quasars: emission lines 



1. INTRODUCTION 



It is widely accepted that the extraordinary activity associated with quasars involves accretion onto a 
supermassive black hole (SMBH). Furthermore, the evidence that almost all massive galaxies host SMBHs 
has become quite convincing. It has been found that SMBH mass is correlated with the host galaxy's bulge 
luminosity (e.g., Kormendy & Richstone 1995; Magorrian et al. 1998; McLure & Dunlop 2001; Marconi & 
Hunt 2003) as well as the stellar velocity dispersion (e.g., Gebhardt et al. 2000a; Merritt & Ferrarese 2001; 
McLure & Dunlop 2002; Tremaine et al. 2002). Because luminous quasars have been observed to reside in 
massive early-type galaxies (McLure et al. 1999; Kukula et al. 2001; McLeod & McLeod 2001; Nolan et al. 
2001; Percival et al. 2001; Dunlop et al. 2003), this implies that the evolution of spheroidal galaxies and 
quasars is intricately tied together (e.g.. Silk & Rees 1998; Haehnelt & Kauffmann 2000; Adams et al. 2001; 
Merritt & Poon 2004; Di Matteo et al. 2005). Therefore, understanding the cosmic evolution of SMBHs is 
an important task of modern astronomy. 
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Reverberation mapping (Blandford & McKee 1982; Peterson 1993) is often used to estimate SMBH 
mass, Mbh, in Type 1 active galaxies (Wandel et al. 1999; Kaspi et al. 2000; Peterson et al. 2004). One 
of the principal advantages of this method is that it does not require high spatial resolution, but rather 
relies on the time lag between the continuum and emission line variability. Under the assumption that the 
broad line region (BLR) is in Keplcrian motion, the time lag is combined with the line width to give an 
estimate of SMBH mass. This, in principal, makes it immediately applicable for both faint and distant 
quasars. Although there are many potential systematic uncertainties in the technique (Krolik 2001), there 
has been good agreement between SMBH masses inferred from reverberation mapping and those inferred 
from the stellar velocity dispersion (Gebhardt et al. 2000b; Ferrarese et al. 2001; Nelson et al. 2004; Onken 
et al. 2004). 

Unfortunately, reverberation mapping requires long-term intensive monitoring, which is not practical 
for large samples. In addition, the long time scales for variability in bright high redshift sources make 
reverberation mapping for these sources unfeasable. Fortunately, a correlation has been found that relates 
the BLR size and monochromatic continuum luminosity (the R-L relationship, Kaspi et al. 2000, 2005), 
making it possible to estimate virial masses from single-epoch spectra using the H/3 width (Wandel et al. 
1999). In addition, Wu et al. (2004) find a relationship between BLR size and the luminosity of H/3, but 
Woo & Urry (2002) do not find any correlation between SMBH mass and bolometric luminosity. H/3 is 
redshiftcd into the near-infrared at z ^ 0.9, making it difficult to observe from the ground for large samples. 
Vestergaard (2002) and McLure & Jarvis (2002) have argued for the use of C IV and Mg II, respectively, to 
estimate virial masses from single-epoch spectra, allowing quasar SMBH masses to be estimated at high z 
from the ground. Many studies have exploited these results and estimated SMBH masses for large samples 
of quasars (e.g., Bechtold et al. 2003; Corbett et al. 2003; Warner et al. 2003; Vestergaard 2004; McLure & 
Dunlop 2004; McLure & Jarvis 2004). Dietrich & Hamann (2004) have found that the estimates based on 
H/3 and C IV agree well for high z quasars, whereas Mg Il-based estimates are typically a factor of 5 times 
lower. In constrast to this, Shemmer et al. (2004) and Baskin & Laor (2005) argue that C IV does not give 
as accurate of an estimate of SMBH mass as H/3. 

Previous methods that utilize single-epoch spectra have relied on empirical linear relationships that 
estimate the BLR size for a given source luminosity, and then use this estimate of R in the virial relationship. 
Although this will give good estimates for R, this may give less efficient estimates for SMBH mass. The 

reason for this is that the luminosity may also be correlated with Mbh in a manner independent of the R-L 
relationship, such as through the accretion process. By only using the luminosity as a proxy for R when 
estimating Mbh, one is ignoring the additional information of Mbh that is contained within L. In other 
words, the standard broad line estimates are based on the distribution of velocities at a given luminosity and 
Mbh, as folded through the R L correlation and imder the assumption of Keplcrian motion. However, the 
distribution of black hole masses at a given luminosity and line width is not simply obtained by inverting 
the distribution of velocities at a given luminosity and black hole mass, unless L and Mbh are statistically 
independent. 

Recent studies have found evidence for a correlation between luminosity and black hole mass (e.g., 
Corbett et al. 2003; Netzer 2003; Peterson et al. 2004), suggesting that black hole mass and luminosity are 
not statistically independent. However, in contrast to other studies. Woo & Urry (2002) have argued that 
there is no significant correlation between bolometric luminosity and Mbh- If a correlation between Mbh 
and L exists, then we can combine the Mbh-L correlation with the broad line mass estimates to obtain, on 
average, more accurate estimates of Mbh for a given luminosity and line width. 

In this work, we outline a formalism that allows one to estimate the probability distribution of an AGN's 
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black hole mass, given some set of observables. We focus on the special case of estimating Mbh given some 
monochromatic luminosity, and the width of an emission line. We also search for other parameters of 
the C IV line that may contribute additional information of Mbh, thus decreasing the uncertainty in the 
estimated Mbh- Although it is possible to include other quasar properties, such as radio loudness, X-ray 
loudness, or variability (e.g., Xic ct al. 2005; O'Neill ct al. 2005; Pcssah 2006), for simplicity wc have chosen 
to only include parameters that may be measured from a single spectrum containing C IV. This allows Mbh 
to be estimated using only one spectrum, and it is thus not necessary to compile observations from several 
different spectral regions. 

We use the formalism developed here to justify using the C IV line width in estimating Mbh, and 
attempt to indirectly infer the C IV R-L relationship. We have chosen the C IV line because it is readily 
observable from the ground over a wide range in redshift (1.5 < 2 < 4.5), is less effected by blends with 
iron and other lines, has shown to give consistent mass estimates with H/3 (Vestergaard 2002), is commonly 
employed to estimate SMBH mass, and archival UV spectra are available for most of the sources with 
reverberation-based masses (Peterson et al. 2004). Despite its common usage, the R-L relationship for C 
IV is mostly unexplored, as there are only seven data points with reliable C IV reverberation mapping data 
(Peterson et al. 2005). Often it is assumed that C IV BLR size has the same dependence on luminosity as 
H/3 (e.g., Vestergaard 2002; Netzcr 2003; Vestergaard & Peterson 2006). Another possibiHty is to assume 
R oc L^/^, as predicted from simple photoionization theory or if the BLR size is set by the dust sublimation 
radius (Netzer & Laor 1993); this was done by Wandel et al. (1999) and Shields et al. (2003) for the H/3 
line. Peterson et al. (2005) performed a linear regression using five AGN with a total of seven data points 
and find R oc L^y similar to that for the Balmer lines^. Unfortunately, this result is almost entirely 

dependent on the inclusion of NGC 4395, the least luminous known Seyfert 1 galaxy, since the data points 
for the other AGN are clustered around ALa(1350A) xq^^t^ ergs s~^. Vestergaard & Peterson (2006) 
found that using the C IV line to calculate black hole masses assuming R oc gave results consistent 
with masses obtained by reverberation mapping. 

The layout of the paper is as follows. In § 2 we describe the general problem of estimating SMBH 
mass from single-epoch spectra, and in § 3 we describe our two samples. The first sample is a set of 25 
quasars with black holes from reverberation mapping and archival UV spectra, and the second sample is a 
set of 100 quasars for which we have spectra containing both the H/J and C IV emission lines. In § 4 we 
describe the method we employ to estimate the emission line parameters. In § 5 we test if any other UV or 
C IV parameters contribute useful information about Mbh- In § 6.1 we use the sample with reverberation 
mapping data to investigate the Mbh-L relationship, and in § 6.2 we test if including the C IV line FWHM 
is preferred by this sample and investigate the nature of the C IV R-L relationship. In § 7 we use our larger 
sample to compare estimates of Mbh obtained from both single-epoch H/3 and C IV. We summarize our 
results in § 8. 

In this work we adopt the WMAP best-fit cosmological parameters {h = 0.71, = 0.27, = 0.73, 
Spergel et al. 2003). We will use the common statistical notation where a point estimate of a quantity is 
denoted by placing a "above it, e.g., Mbh would be an estimate of Mbh- We will also commonly refer to the 
bias of an estimate, where the bias of Mbh is Bias = E{Mbh) ~ Mbh- Here, E(Mbh) is the expectation 
value of Mbh- An unbiased estimate of Mbh is one such that E{Mbh) = Mbh- 



^This value appears in an erratum to this paper 
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2. ESTIMATING BLACK HOLE VIRIAL MASSES 

If one assumes that the BLR gas is in Keplerian motion, then the mass of the central black hole, Mbh, 
may be estimated from the virial theorem : 

MBH = f-^. (1) 

Here R is the distance between the BLR that is emitting a particular line and the central continuum source, 
V is the velocity dispersion of the line-emitting gas, and G is the gravitational constant. The velocity 
dispersion is inferred from the width of the line, quantified using either the FWHM or line dispersion, cr*. 
The factor / is a scaling factor that converts the virial product, RV'^ /G, to a mass. Typically, / has been set 

to the value appropriate for an isotropic velocity field, where / = 3/4 if one uses the broad line FWHM to 
estimate V (e.g., Netzer 1990); however, there may be systematic effects that can significantly effect the value 
of / (Krolik 2001). Onken et al. (2004) used the correlation between Mbh and stellar velocity dispersion to 
estimate an average scale factor of (/) = 5.5 when the line dispersion is used to estimate V . For simplicity, 
in the rest of this work we will assume / = 1, so that what is really being estimated is the virial product. 
After estimating the virial product, we can convert it to a mass using any adopted value of /. 



2.1. Estimating Mbh from Reverberation Mapping 

For the case of reverberation mapping, estimating Mbh is straightforward. The BLR size, R, can be 

estimated as ct, where c is the speed of light and r is cither the peak or the centroid of the line-continuum 
cross-correlation function. The velocity dispersion of the line-emitting gas, V, is estimated from the width 
of the broad emission line as measured in the variable part of the spectrum. Measuring V from the variable 
(RMS) spectrum ensures that one is probing the line emission that is actually varying, i.e., the BLR gas 
that is at the distance ct. One then uses cr and V in Equation (1) to estimate Mbh- See Peterson et al. 
(2004) for recent reverberation mapping results. 



2.2. Estimating Mbh from Single-Epoch Spectra 

Estimating Mbh from a single-epoch spectrum (SES) is a somewhat different problem for several rea- 
sons. In this section we illustrate why the SES case is different, and provide a general formalism for calculating 
the conditional probability distribution of Mbh, given a set of obscrvablcs. 

In the single-epoch case, one cannot directly observed the BLR size, but instead employs a correlation 
between R and continuum luminosity. However, there is considerable scatter in the R-L relationship, which 
is propogated through when using L instead of R to estimate Mbh- In this case, the conditional probability 

distribution of Mbh- given the broad line estimate, is broad, typically a factor of a few (Vestergaard 2002; 
Kaspi et al. 2005). This scatter can be reduced by incorporating information about Mbh from other 
observables. In the more general case, one may have several parameters that contain information about 
Mbh, such as emission line width, variability, luminosity, etc. For example, Merloni et al. (2003) have 
found evidence that black hole mass is correlated with X-ray and radio luminosity. In this case, one could 
combine the information from the line width, UV luminosity. X-ray luminosity, and radio luminosity to 
obtain an estimate of Mbh that is more accurate than would have be obtained solely from some subset of 
these parameters. 
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To be more specific, suppose that one has a set of observables, denoted by X. Then, given the observ- 
ables, X, the conditional probabihty of Mbh given X is given by Bayes' Formula 

p(M,.|X) = ^^-^'^;g(^-^) , (2) 

where p{Mbh) is the intrinsic probability distribution of Mbh for a sample, and p{X) is the distribution of 
the observables; p{X) is just a normalizing constant and may be ignored. 

In this work we are concerned with the distribution of Mbh given L, the emission line width, and 
redshift. For this case. Equation (2) becomes 

p{m\l, V, z) oc p{v\l, m, z)p{l\m, z)p{m\z). (3) 

Here, we are using the notation, m = log Mbh, I = logAL;^, and v = logV". The first term, p{v\l,m, z), 
is the distribution of line widths at a given I, m, and z. The second term, p{l\m, z), is the distribution of 
luminosities at a given m and z. The last term, p{m\z), is the distribution of m at a given redshift. As 
we will show later in this section, p{v\l,m,z) is obtained by plugging the R-L relationship into the Virial 
theorem (Eq.[l]). This is the standard method of estimating Mbh from the broad lines, but it implicitly 
assumes p{m\l, v, z) oc p{v\l, m, z), and therefore that p{m\z) and p{l\m, z) are uniform. Taking p{m\z) and 
p{l\m, z) to be uniform results in a broader distribution of Mbh, given L, V, and z, and thus a less efficient, 
albeit still unbiased, estimate of Mbh- However, by incorporating information on both the distribution of 
luminosities at a given black hole mass and redshift, and the distribution of Mbh at a given redshift, one 
can obtain a better estimate of Mbh- 

To estimate a functional form for p{l\m) and p{v\l,m, z), suppose we observe some quasar with SMBH 

mass Mbh at a redshift z, where m is drawn from some probability density p(m|z), m|z ~ p(rn\z). We 
assume that the accretion process for this source generates a luminosity, Lx oc M^]^^^^ , by 

l\m,z = aQ + am{z)m + ei{z,m). (4) 

Hero, ao is some constant of proportionality, and ei{z.m) is the random error term representing the scatter 
about this relationship. We will refer to Equation (4) as the Mbh-L relationship. The stochastic term, 
ei{z,m), encompasses variations at a given Mbh in the accretion rate (M), radiative efficiency (e), source 
inclination, bolometric correction {Cboi), etc.; since we do not observe these quantities we model them as 
being random. For now, we allow the logarithmic slope, am, to depend on z, and the random error to depend 
on z and Mbh- 

The parameter can be predicted from accretion physics. The radiated bolometric luminosity from 
an accretion fiow can be written as 

L = (1.26 X 10^')e"i^^ W, (5) 

where m = M/Medd is the accretion rate normalized to Eddington. If one assumes Lx = CboiL, then from 
Equation (5) it follows that for this case am = 1 and ei ~ log e + log to + log Ctoi • A more careful analysis of 
the thin disk case suggests XLx oc {MbhM)°'"^ ,am = 2/3, after employing some simplifying approximations 
(e.g., Collin et al. 2002). 

Given this luminosity, the BLR distance R is assumed to be set by the luminosity according to the R-L 
relationship, R oc i^' : 

r\l,z = eo + eil + er{z). (6) 
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Similar to above, r = \ogR and er{z) is the stochastic component. In Equation (6) we assume that given 
L, R is independent of Mbh- In addition, we have not assumed any redshift dependence for 9i because it is 
likely that the form of the R-L relationship is independent of z (Vestergaard 2004). However, in Equation 
(6) we have allowed for the possibility of a redshift dependence for as the the intrinsic scatter in the R-L 
relationship may depend on z. 

1 /2 

Simple photoionization theory predicts that R oc L^j^^^, where Lion is the luminosity of the ionizing 
continuum (e.g., Wandel et al. 1999; Kaspi et al. 2000). From the definition of the ionization parameter, U, 
it follows that 

2r = log Lion - log(47rc) -\ogU - log rig - log E, (7) 

where E is the average energy of an ionizing photon and rig is the BLR gas density. If we make the simplifying 

assumptions that r is set by Equation (7), L\ oc Liom and the means of the distributions of /7, rig, and E are 
independent of La, then comparison with Equation (6) shows that 9i = 1/2 and = — log C/ — log rig — log .B. 
If this is not the case, but rather the means of the distributions of U, Ue, and E have a power-law dependence 
on Lion, then Equation (7) is still valid, but in general 6i ^ 1/2. Either way, in this model the scatter about 
the R-L relationship is partly the result of variations in U, n,,, and E. Other sources of scatter may include 
variations in the conversion between L\ and Lion, source inclination, and the non- instantaneous response of 
the BLR to continuum variations. 

If R is set by the dust sublimation radius, then we also expect 9i = 1/2 (Netzer & Laor 1993), but in 
this case R oc L^/^. A relationship of the form R oa L^J^"^ is consistent with the results of Peterson et al. 
(2004) for the Balmer lines if one uses the UV continuum luminosity. 

Finally, from Equation (1) the observed SES line width depends on R and Mbh as 

1 1 

v\r,m = VQ- ^ (8) 

where, = \o<g{\fG j J ses) , and is the stochastic term. The term fsES converts the SES line width 
measurement into a velocity dispersion, and e„ describes random deviations of single-epoch v from that for 
the RMS spectrum. Vestergaard (2002) has found that the single-epoch H/3 FWHM and the RMS H/3 
FWHM arc consistent with a scatter of ^ 20%, if one does not subtract the SES H/? narrow component. 
When the FWHM is used, fsES ~ 1/2. Through the stochastic terms ei,er, and Cy, Equations (4), (6), 
and (8) define the conditional probability densities that describe how L and V depend on Mbh, and how R 
depends on L. 

It is useful to examine the special case of Gaussian error terms, €i,er, and and Gaussian p{m\z). 
In this case, the distribution of m at a given I and v is also normal. If we assume that ei,er, and e„ are 
uncorrelated, have zero mean, and variances (Tf,<T^, and cr^, respectively, then the optimal estimate of m 
is the mean oi p{m\l,v,z), ji. If we make the further assumption that ei,er, and Cy are independent of m 
and z, and that does not depend on z, then p{m\l,v,z) takes a particularly simple form. The mean of 
p{m\l, V, z) may be calculated from the properties of the normal distribution (e.g., Gelman et al. 2004) as 

^ mBL/<J%L + mML/crliL + TO(z)/cr^(2:) 
^ l/a%^ + l/al,^ + l/a^z) " 

Here, niBL = m^^ + 9il + 2v is the standard broad line mass estimate, found by plugging Equation (6) into 
Equation (8), = ^ "^vq, Ubl = + 4ct^ is the variance in ttibl, itiml = {I — cto)/ctm is the estimate 
of m based on the Mbh-L relationship, a^j^ = af /a^ is the variance in niML, fh{z) is the mean m at a 
given z, and <t'^{z) is the variance in m at a given z. 
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The variance in n as an estimate for m, ct^, is 
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(10) 



As is apparent from Equation (10), the variance in /x as an estimate for m is always less than that for the 
standard broad line estimate, tjibl- In fact, in the limit (afjj^,a'^{z)) oo, p(?7?,|Z) and p{m,\z) supply no 
information on the black hole mass and jj, converges to the broad line estimate. However, by combining the 
broad line mass estimate with the information on Mbh that is contained within the luminosity and redshift 
we can obtained a better estimate of m. 

In reality, the distribution of m is unlikely to be Gaussian, and since m is a component of e;, this 

would violate the assumption that is Gaussian. Instead, the distribution of m is likely bimodal (Ho 2002; 
Marchesini et al. 2004; Hopkins et al. 2006; Cao & Xu 2006) as a result of a transition from a radiatively 
inefficient flow to an efficient one at m ~ 0.01 (e.g.. Jester 2005). However, because of the flux limits of 
modern surveys, most observed broad line quasars will have m > 0.01 (McLure & Dunlop 2004; Vestergaard 
2004). This results in a imimodal, relatively smooth and symmetric distribution of m for observed quasars 
(Hopkins et al. 2006). It may also be that BLRs do not form in sources with rh < 10~^ (Nicastro et al. 
2003; Gzerny et al. 2004), and therefore the distribution of rh for broad line AGN would have a lower limit 
at m 0.001. If true, then unimodality in the distribution of rh for broad line AGN is ensured. 

The stochastic scatter about the Mbh-L relationship, e^, is the sum of random deviations in rh, bolo- 
metric correction, radiative efficiency, etc. We assume that the distributions of most, if not all, of the 
constituent components of ti are not too different from a Gaussian distribution, i.e., unimodal, smooth, and 
fairly symmetric. Then, by the central limit theorem, the distribution of the sum of these components, ei, 
will tend toward a normal density. Therefore, without any evidence to the contrary, the normal density 
should provide an accurate approximation to the true form oip{l\m). A similar argument may be used to 
justify the assumption of normality for and e^. 

To illustrate the improvement that Equation (9) offers over the broad line mass estimate, we simulate 
values of m, I, and v. The simulations were performed as follows. First, we simulated values of m from a 
smoothly-connected double power-law, with a mean of m « 7.87 and a dispersion of cr « 0.45 dcx. Values 
of / were then generated according to Equation (4), with — 37, a^n = 1- The gaussian scatter about this 
relationship had a dispersion of 0.7 dex. Then, we generated values of v as v = 11 — //4 -|- m/2 -|- e, where e 
was a gaussian random variable with dispersion 0.2 dex. This form for v assumes R oc V-/"^, and corresponds 
to a intrinsic scatter in the broad line mass estimates of 0.4 dex. The parameters for simulating / and v were 
chosen to be similar to the results found in § 6. Finally, we calculated broad line mass estimates from the 
simulated luminosities and line widths, rhsL = —22 + 1/2 + 2v, and mass estimates fi according to Equation 



The results are shown in Figure 1. As can be seen from the distribution of residuals, the mass estimates 
that combine all available information on the black hole mass, /i, are more accurate on average than the 
broad line estimates, rhsL- In addition, the distribution of n provides a more accurate estimate of the true 
distribution of m than does ttibl, with the distribution of m inferred from rfiBL being too broad. It is 
interesting to note that both of these results are in spite of the fact that the intrinsic distribution of m is not 
Gaussian, which was assumed when deriving Equation (9). This suggests that if the intrinsic distribution of 
m for a sample is not too different from a normal density. Equation (9) will still give more efficient estimates 
than the broad line estimates. This is reasonable, considering that Equation (9) 'shrinks' the black hole 
mass estimates towards the sample mean by an amount inversely proportional to the intrinsic variance in m 



(9). 



- 8 - 



of a sample. 



2.3. Cautions for Using Quantities Calculated from the SES Mass Estimates 

The intrinsic uncertainty in Mbh inferred from Lx and the hne width may be thought of as the 'mea- 
surement' error in Mbh- This intrinsic uncertainty can cause problems when using the black hole mass 
estimates, m, to calculate additional quantities. In particular, quantities based on the square of the esti- 
mated m, such as correlation coefficients and linear regressions, can be significantly effected (e.g., see Akritas 
& Bcrshady 1996; Fox 1997). 

Suppose we are interested in calculating the correlation between Mbh and some other parameter X. To 
estimate this correlation, we would obtain a sample of quasars with black hole masses estimated from some 
assumed form oi p{m\l,v, z), rh. Typically, these are the standard broad line estimates. We then calculate 
the correlation coefficient between X and rh. However, because we do not have the actual black hole masses 
for our sample, but instead obtained estimates from the continuum luminosities and widths of the broad 
lines, this is not the true correlation coefficient between X and m. Because our estimated black hole masses 
have been 'measured' with error, this broadens the observed distribution of black hole masses, and thus 
biases the observed correlation coefficient towards zero. Therefore, the correlation coefficient obtained from 
the estimated black hole masses will be, on average, less in magnitude than the correlation coefficient that 
would have been obtained using the actual black hole masses. 

To prove this point, we note that the logarithmic black hole mass estimates are related to the actual black 

hole masses as m = m + e„j, where is the random error between m and m. The linear correlation between 
the parameter of interest X and m is p = Covar{X,m)/[Var{X)Var{m)]^^^ , where Covar and Var are the 
sample covariance and variance, respectively. Since we don't actually observe to, we can't calculate the true 
correlation coefficent, but instead we calculate the correlation between X and m, p. While the covariance 
between X and to is unaffected by using to instead of to, the sample variance of to is Var{m) = Var{rn)-\-a'^, 
where cr?^ is the intrinsic uncertainty in rh as an estimate for to, ct^ = E{e^). The observed correlation is 
then given by 

^^^(^) (11) 

Var{m) + al-\ ^ 

Therefore, correlation coefficients calculated using the estimated black hole masses will be reduced in mag- 
nitude from the true correlation by a factor of [l -|- a^/Var{m.)] . For broad line mass estimates based 
on H/3, Vestergaard & Peterson (2006) find = 0.43 dex. In this work we find that cta ~ 0.40 dex for 
C IV-based broad line estimates . For a sample with an intrinsic dispersion in to of 0.75 dex, these values 
of (T,jj correspond to a decrease in the magnitude of any observed correlation by « 12%. If the sample has 
an intrinsic dispersion of 0.4 dex, similar to the intrinsic uncertainties in the broad line estimates, then the 
magnitude of the observed correlation coefficient is reduced by w 30%. These effects is not negligible, and 
can be more serious for linear regression (e.g.. Fox 1997). In light of these issues, care must be taken when 
calculating quantities from black hole mass estimates based on single-epoch spectra. 

Another problem arises when one is using a flux limited sample to estimate the intrinsic distribution 
of the black hole mass, i.e., the active black hole mass function (Wang et al. 2006), based on broad line 
estimates. There is currently significant interest in this problem, as the active black hole mass function is 
an important tool in understanding SMBH formation and civolution. Unfortunately, the limiting flux of a 
survey causes incompleteness in black hole mass, the degree of which depends on the distribution of i at a 
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given TO and z, p(l\m,z). In addition, because the broad line mass estimates arc measured with error, the 
distribution of m inferred from the broad hnes can be significantly broader than the intrinsic distribution 
(cf., Fig.l). Because of these issues, estimates of the active black hole mass function obtained from broad 
line estimates should be interpreted with caution. 

In order to estimate the degree of incompleteness in m, it is necessary to re-express the survey selection 
function as a function of to. Following Gelman et al. (2004), we introduce an indicator variable, /, denoting 
whether a source is included in the survey, where /= 1 if a is source included and / = if a source is missed. 
Then the selection function of the survey is the probability that a source is included in the survey for a given 
luminosity and redshift, p{I = z). The selection function for black hole mass is then 



As can be seen, estimating the completeness in m for a survey depends on the form of p{l\m, z). Therefore, 
it is important to understand p{l\m, z), even if such an understanding does not result in significantly better 
estimates of to. 



We employ two samples in our analysis. The first sample consists of a set of 25 low-z sources that have 
reverberation mapping data from Peterson et al. (2004). Wc use this sample to investigate the Mbh^L 
relationship and the C IV R-L relationship, and to justify using the C IV line to estimate AGN black hole 
masses. The second sample is a set of 100 quasars for which we have seperate spectra containing the H/J 
and C IV emission lines. We use this sample to compare the broad line estimates of Mbh obtained from the 
two lines. 

Peterson et al. (2004) calculated virial products for 35 AGNs based on the reverberation mapping 
method. Of those 35 sources, we selected ones with archival UV spectra. We did not include those sources 
which Peterson et al. (2004) listed as having unrehable virial products (PG 0844+349, PG 1211+143, PG 
1229+204, and NGC 4593). We also did not include NGC 3227 or NGC 4151, as the C IV hne for these 
sources had significant absorption. In addition, we removed NGC 4051 from the analysis because this source 
is an outlier in the BLR R-L relationship (Vestergaard 2002; Kaspi et al. 2005), and thus probably does not 
follow the linear relationship assumed in Equation (6). Our sample consists of 25 AGN: 14 sources with HST 
FOS spectra, 9 sources with lUE spectra, 1 source with HST GHRS spectra, and 1 source with HST STIS 
spectra. Of the 14 sources with FOS spectra, 13 were taken from Bechtold et al. (2002), and the other (NGC 
5548) from Evans & Koratkar (2004). For archival sources with more than one observation, we took the 
source with the longest exposure time. All spectra are single-epoch, except for NGC 5548, which is averaged 
over approximately a month of HST observations. The sample is summarized in Table 1. 

The luminosities were calculated from the predicted continuum flux at 1549A, assuming a power law 
continuum (see § 4.1). We corrected luminosities for galactic absorption using the E{B — V) values taken 
from Schlegel et al. (1998), as Usted in the NASA/IPAC Extragalactic Database (NED), and the extinction 
curve of Cardelli et al. (1989), assuming a value of Av/E{B — V)= 3.1. We did not do this for the Bechtold 
et al. (2002) sources as they have already been corrected for galactic absorption. The C IV line widths 
for the lUE sources were corrected by subtracting an assumed instrumental resolution of 1000 km s^^ in 
quadrature from the measured line widths; the resolution for the other instruments is negligible compared 
to the emission line widths, so no correction was performed. 




(12) 



3. 



THE SAMPLE 
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Wc compiled UV and optical spectra for a sample of 100 sources for the purpose of comparing the C 
IV-based estimates derived here with the estimates based on H/3 and the empirical R-L relationship. Of 
these sources, 89 have z < 0.8, 6 have z 2.3, and 2 have z ~ 3.3. The UV spectra for the z < 0.8 sources 
are FOS spectra from Bechtold et al. (2002), the optical spectra for 51 of these sources are from Marziani 
ct al. (2003), and the optical spectra for 9 of these sources arc from the SDSS DR2 (Abazajian et al. 2004). 
Optical spectra for the 29 remaining z < 0.8 sources were obtained by us. Twenty-seven of the sources were 
observed at the Steward Observatory 2.3m Bok Telescope on Kitt Peak using the 600 lines mm~^ grating 
of the B&C Spectrograph; these sources had moderate spectral resolution (~ 5 A). The other two sources 
were observed at the Magellan Baade Telescope using the Inamori Magellan Areal Camera and Spectrograph 
(IMACS); long-slit spectra were obtained for these sources using a slit width of 0.9" in long camera mode. 
We used the 600 lines mm"! grating for PKS 1451-375 and the 300 lines mm"! grating for PKS 2352-342, 
giving spectral resolutions of 2 A and ~ 5 A, respectively. The log of new spectra is displayed in Table 2, 
and they are shown in Figure 4. The spectra were reduced using the standard IRAF routines. 

Rest-frame UV spectra for the eight sources with z > 2 arc from Scott et al. (2000) and the SDSS. The 
continuum luminosities and H/3 FWHM for the z ~ 2.3 sources were taken from Mcintosh et al. (1999) 
and corrected to our adopted cosmology. The continuum fluxes for the two 2; 3.3 sources were measured 
directly off the published spectra of Dietrich et al. (2002), and the values of H/3 FWHM for them are from 
Dietrich k Hamann (2004). 

4. LINE PROFILE PARAMETERS 

4.1. Extracting The Line Profile 

In order to extract the C IV emission line, it is necessary to subtract the continuum, Fe emission, and 
the He II A1640 and O HI] A1665 emission lines. To remove the continuum and iron emission, we use a 
variation of the method outlined in Boroson & Green (1992). We model the continuum as a power law 

of the form /,y oc v°' . The Fc emission was modeled as a scaled and broadened form of the Fc emission 
template extracted from I Zw I by Vestergaard & Wilkes (2001). In constrast to most previous studies, we 
simultaneously fit the continuum and Fe emission parameters using the Levenberg-Marquardt method for 
nonlinear x^-minimization; this is similar to the method used by Mcintosh et al. (1999). Performing the fits 
in this manner has the advantage of providing an estimate of the uncertainties in these parameters, given by 
the inverse of the curvature matrix of the space. The set of possible windows used to fit the continuum 
and Fe emission are shown in Table 3. The actual continuum windows used to fit the continuum and Fe 
emission for any particular source depended on that source's available spectral range. 

Having obtained an estimate of the continuum and Fe emission, wc subracted these components from 
each spectrum. Wc then extracted the region within —2 x 10^ km s~^ and 3 x lO** km s^^ of 1549A. Here, and 
throughout this work, we will use the convention that negative velocities are blueward of a given wavelength. 
Narrow absorption lines were removed and interpolated over. We then removed the He II A1640 and O HI] 
A1665 emission lines from the wings of the C IV profile. This was done by modelling the C IV, He II, and 
O III] lines as a sum of Gaussians. In general, C IV was modelled as a sum of three Gaussians, He II two 
Gaussians, and O III] a single Gaussian, however this varied from source to source. The C IV extraction 
was done interactively in order to ensure accuracy of the fits. After obtaining estimates of the He II and O 
HI] profiles, we subtracted these components. We did not fit the N IV] A1486 emission line as this line is 
typically weak and lost in the C IV wings. 



Extraction of the H/? profile was done in a similar manner. For the optical Fe emission we used the I 
Zw I template of Veron-Cetty et al. (2004). After subtracting the continuum and Fe emission, we extracted 
the region within ±2 x 10^ km s"^ of 486ll. The H/3 profile was modeled as a sum of 2-3 Gaussians. The 
[O III] A4959A and [O III] A5007A lines were modeled as a sum of 1-2 Gaussians, depending on the signal- 
to- noise of the lines. A sum of two Gaussians was used for the higher S/N lines because the [O III] lines 
are not exactly Gaussian, and not because the seperate Gaussians are considered to be physically distinct 
components. The widths of the narrow Gaussians for [O III] were fixed to be equal. The [O III] lines and the 
narrow II/3 line were then subtracted to extract the II/3 line. As with the C IV profile, the entire extraction 
was done interactively. 



4.2. Estimating The Line Profile Parameters 

We measured the C IV line shift, FWHM, EW, and the first and second line moments. We define 
the line shift, Av, as the location of the line peak relative to 1549A, in km s~^. The location of 1549A is 

determined from the redshifts listed on NED. We have checked the references for the redshifts to investigate 
how z was estimated, but not all of the references report this. For those that did specify, the redshifts were 
determined from the narrow emission lines (e.g., [O III] A5007.). 

lUE observations were done using a very large aperture, so the values of Av for the lUE sources may be 
biased. However, there is no noticeable difference between Av estimated from the lUE spectra, and those 
estimated from the HST spectra. Furthermore, the Av parameter docs not enter into our final analysis, and 
so even if the lUE Av parameters are significantly biased our conclusions remain unchanged. 

The first moment of the C IV line is the centroid, nciv, and the square root of the second central 
moment is the line dispersion, a* . The zeroth line moment is the line flux. The line moments are calculated 
as 



F = 


n 
i=l 




(13) 


l^civ = 


n 

i=i / 


i=l 


(14) 


crl = 


n 

i=i / 


1=1 


(15) 



Here, F is the line flux, (5A is the spacing between subsequent wavelengths, y is the observed spectral flux 
density, and x is the velocity relative to 1549A. In practice we do not perform the sums over all the n 
data points, but rather only over those data points with /j > 0.05(max/), and assuming that the profile 
is monotonically decreasing blueward and redward of the peak. Here, / denotes the best-fit line proflle, 
found from modelling the emission lines as a sum of Gaussians. This allows us to define the extent of the 
line profile. Although this gives us biased measurements of the line moments, it keeps these measurements 
stable. In particular, a* can be very sensitive to the profile wings, which have the highest uncertainty. By 
truncating the line proflle we keep the estimate stable and less sensitive to errors in the Fe and continuum 
subtraction, as well as errors in the He II A1640 and O III] A1665 subtraction. 

The line moments are calculated using the observed line profile, y, and not the best-flt to the spectral 
flux densities, /. We do this because the line moments are relatively insensitive to a lack of smoothness. 
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To be specific, consider a line profile, f{x), and its Fourier transform, /(fc). For simplicity, we consider 
the continuous case here. The j*'' unnormalized line moment, may be written in terms of the Fourier 
transform of the line profile: 

,,=y_.^/(.)dx=^, (16) 

where f^-'\k) is the j*-^ derivative of f{k) and i ~ \/— 1. One can sec from Equation (16) that the line 
moments only depend on f{k) near fc = 0, and are thus insensitive to the high frequency behavior of f{x). A 
generic smoothing operator, such as a Gaussian fit, will shrink the high k components of /(fc) more than the 
low k, as it is generally the case that the high k components have lower signal-to-noise. However, because the 
line moments do not depend on the high k components, nothing is gained by enforcing smoothness. Because 
of this we just use the observed C IV profile, y, as it is an unbiased estimate of the true profile, /, whereas 
the best-fit estimate, /, is a biased estimate. 

There has been some discussion in the literature over whether FWHM or ct* is the better width to use 
in estimating m (Fromerth & Melia 2000; Peterson et al. 2004). Promerth & Melia (2000) suggested using 

the line dispersion, cr, , arguing that it provides a better estimate of the velocity dispersion for an arbitrary 
line profile. Other advantages of u* include its insensitivity to noise and narrow absorption lines. However, 
cr* can be significantly biased due to its sensitivity to blending with other lines in the line wings, and to the 
removal of continuum and iron emission. These facts can be understood in light of Equation (16), which 
shows that the second moment of the line depends on the low-fc behavior. Because of this, (T* is relatively 
unaffected by information on small scales, such as noise and narrow absorption lines, but is significantly 
affected by information on large scales, such as line blending, truncation of the line profile, and continuum 
placement. The FWHM, on the other hand, is insensitive to line blends and errors in the continuum and 
iron subtraction. In addition, FW HM is easily measured, however it probably provides a poor estimate of 
the velocity dispersion for irregularly shaped lines profiles. Peterson et al. (2004) compared the strengths 
and weaknesses of these two measurements and concluded that cr* was the better parameter when measured 
in the RMS spectrum. 

In order to choose the better SES estimate of the BLR velocity dispersion, we calculate the partial 
correlation between m and FWHM and cr*, respectively. The partial correlation coefficient describes the 
correlation between m and line width at a given luminosity; the line width with the higher partial correlation 
should give a better estimate of m. The partial linear correlation between log FWHM and m is 0.36, while 
the partial linear correlation between loga* and m, is 0.31. Because the FWHM has a moderately higher 
partial correlation, and because the FWHM is not as affected by errors in the line deblending, continuum 
placement, etc., we use the FWHM as an estimate of the velocity dispersion throughout the rest of this 
analysis. 

We only measure the FWHM of the H/3 emission line. This is because we are only concerned with 

getting an estimate of Mbh from single-epoch H/3 based on the virial theorem for comparison with our C 
IV-based estimates, so the only line parameter of interest is the H/3 width. 

The standard errors on the line parameter measurments are estimated using the bootstrap (Efron 1979). 
In this method, we take our best-fit line profile spectral fiux densities, /, and generate riboot = 128 simulated 
observed line profiles by adding Gaussian noise to /. We then estimate the line parameters of the simulated 
line profiles and calculate the variance in these parameters over the bootstrap samples. 

In Figure 3 we plot logMBH against logALA(1549 A), logFWHM,logait,logEW, Av, nciv, and con- 
tinuum spectral slope. We report our measurements in Table 4. We have compared our measurements with 
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Wang ct al. (1996), Bechtold et al. (2002), Baskin & Laor (2004), and Kuraszkiewicz et al. (2004), and find 
them to be consistent after accounting for the different procedures used to measure these quantities. For 
clarity, we have removed three outliers in a from the plot of m against a. These sources were 3C 390.3, 
NGC 3516, and NGC 7469. 3C 390.3 is a broad-line radio galaxy with highly variable Balmer lines and 
doublc-pcakcd Ha (Corbett et al. 1998) and H/3 (Ostcrbrock ct al. 1976) profiles; the C IV emission line also 
exhibits a double-peaked profile. The other two sources represent some of the faintest sources in our sample, 
and their spectra may have a contribution from their host galaxies; however, we notice nothing unusual 
about their spectra, save for their unusually soft values of a. 

5. SELECTING THE IMPORTANT PARAMETERS 

It is useful to investigate whether any of the additional parameters that we measured for C IV are 
significantly correlated with m, and thus contribute information about m in addition to that in I and 
FWHM. We used our sample of 25 sources with Mbh from reverberation mapping to test if including 
log EW, a, ficiv, or A?; resulted in more accurate estimates of Mbh- However, we did not find any evidence 
to warrant the inclusion of these parameters. 

In order to assess whether the data support including any additional parameters, we express m as a 
linear combination of all possible combinations of log XLx, log FWHM, log EW, a, jiciv, and Aw, a total of 
2^ = 64 subsets. The regression coefficients are estimated via least-squares, and the result is a set of 64 
linear regressions. In order to assess the relative merits of each of the regressions, and thus each of the 
respective parameters, we employ the Bayesian Information Criterion {BIC Schwartz 1979). Using the BIC 
allows us to undertake a Bayesian comparison of the models without actually carrying out the full Bayesian 
prescription, considerably simplifying things. The BIC is easily obtained from the log-likelihood of the data 
as 

BIC = 2t{i)) + dlnn. (17) 

Here, is the log-likelihood of the data evaluated at the maximum-likelihood estimate, V denotes the 
regression parameters, d is the number of parameters in the regression, and n is the number of data points. 

The only parameter significantly correlated with mis I, with a posterior probability of p; = 0.975. The 
data are ambiguous as to whether FWHM is related to Mbh {pfwhm = 0.562); however, when using the 
BIC we did not compare with regressions that assume Mbh FWHMqjy, and we perform a more careful 
analysis in § 6.2. In addition, the data give weak evidence that the remaining parameters are unrelated to 
m. The posterior probabilities that m is correlated with these parameters are estimated to be 0.270, 0.210, 
0.198, and 0.186 for logEW, nciv, ^v, and a, respectively. 

6. REGRESSION ANALYSIS 

Now that we have ruled out including any additional parameters for estimating Mbh from SES, we 
proceed to investigate the Mbh-L relationship and the C IV R-L relationship. Because p{l\m) and p{v\l, m) 
are statistically independent in their parameters, we can analyze each one seperately. Throughout this section 
we will be using our sample of 25 sources with black hole mass measurements from reverberation mapping. 
We begin by investigating the Mbh-L relationship. 
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6.1. The Mbh-L Relationship 

We fit a linear relationship of the form I = ao + amm, assuming that the scatter about this relationship 

is independent of m and z. Because we arc fitting the distribution of Z at a given m, we use the BCES(F|X) 
(Akritas & Bcrshady 1996) regression. The BCES methods take into account measurement errors in both 
coordinates by correcting their moments; however, the measurement errors are small compared to the variance 
in both m and I, so the correction is small. Based on the regression, we find 

I = 35.72(±1.67) + 1.17(±0.22)m. (18) 

Here, I = log ALa(1549A). The empirical value of ctm = 1.17 ± 0.22 is consistent with L oc Mbh, if one 
assumes that L\ oc L. The intrinsic scatter in Equation (18) is estimated to be a; = 0.61 dex. The residuals 

and their cumulative distribution function (CDF) are shown in Figure 4. A Kolmogorov-Smirnov test found 
that the regression residuals arc not significantly different from a normal distribution, implying that p{l\m) 
is normal with mean given by Equation (18) and standard deviation ai « 0.61 dex. 

The value of dm found here is consistent with Mbh oc iO-9±o i5 found by Netzer (2003) and Mbh oc 
^o.97±o.i6 fo^j^jj Corbett et al. (2003). Peterson et al. (2004) used the 35 AGN with reverberation mapping 
data to estimate the Mbh^L relationship for L\ at 5100A. Using the BCES bisector regression, they find 
Mbh oc L^iqq"^'^^ ■ This is shallower than our result of Mbh oc -^^1549^°'^^, although the two logarithmic 
slopes are consistent within the errors. The difference in the two values most likely results from the different 
regressions used. We used the BCES(Y'|X) because we are modelling p{l\m), where as the bisector slope 
gives the regression that bisects the distribution of Z at a given m, and of m at a given I. 

Woo & Urry (2002) did not find any evidence for a correlation between Mbh and the bolomctric 
luminosity. They used both a sample of reverberation-mapped AGNs and Seyfert galaxies with black hole 
masses derived from the stellar velocity dispersion. In neither case did they find evidence for a correlation, 
in constrast to the results found here and by others. Unfortunately, Woo & Urry (2002) did not perform 
a regression analysis or report a correlation coefficient, so it is difficult to do a quantitative comparison of 
their results with ours. 

We have found here that = 1.17±0.22 for the Peterson et al. (2004) sample. However, these sources 
are all at low redshift and have m = 0.01-1 (Woo & Urry 2002; Peterson et al. 2004), and thus this value of 
am may not be valid for sources with m < 0.01 and z > 0.2. Furthermore, most high z broad line quasars 

have luminosities and AIbh greater than that of the Peterson et al. (2004) sample, and it is possible that 
am 7^ 1 outside of the reverberation mapping sample range. It is also possible that the scatter about the 
Mbh-L relationship depends on z. The most likely source of a redshift dependence for ei is evolution of 
the distribution of m (e.g., Haiman & Menou 2000; Merloni 2004; Hopkins et al. 2006; Steed & Weinberg 
2006), and possibly evolution of the bolomctric correction. The average m has been observed to increase 
with increasing z (McLure & Dunlop 2004), but this is likely due to selection effects. Because of these issues, 
one should exhibit caution when applying this Mbh-L relationship to high z sources. 



6.2. Inferring the C IV R—L Relationship 

In this subsection, we investigate the C IV R-L relationship using the reverberation mapping sample. 
We perform a linear regression to fit for the value of 6';,i? oc by noting that Mbh oc L'^^V'^ . Based on 
the results obtained here, we find Mbh oc L'-^^I^^-^'^ FWEM^jy, consistent with Mbh oc L\'^lgFWHMl;jy 
expected from simple photoionization theory. 
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To estimate 9i, we fit a linear relationship of the form 2v = m — rriQ^ — Oil using the BCES(y|X) 
method. Here, the free parameters are and 9i. The result is 

rhciv = -17.82(±2.99) + 0.41(±0.07)/ + 2\ogFWHMciv- (19) 

The intrinsic scatter about this relationship is aciv = 0.40 dex, where we have corrected for the measurement 
errors in m and FWHMciv- The correlation between the regression coefficients is Corr{rhQ^^ ,9i) = 
—0.9996. The regression results are consistent with the expectation from simple photoionization theory, 

1 /2 

R cx i;^54g (cf. Eq.[7]). A Kolmogorov-Smirnov test found that the regression residuals are not significantly 
different from a normal distribution. 

Performing the regression with v as the dependent variable ensures that the Mbh-L relationship does 
not 'absorb' into the regression coefficients, as the regression models the distribution of f at a given I and 
m. However, if wc had performed the regression by fitting the distribution of m as a function of I and 
V, the intrinsic correlation between m and I would have absorbed into the results. In this case we would 
be modelling p{m\l,v), and the statistical model would be over-parameterized because we have only three 
variables. However, by sepcrating p{m\l^v) into p{l\m) and p{v\l^'m), as was done in § 2.2, we can analyze 
each distribution seperately and uniquely determine their parameters. 

The residuals of the C IV broad line estimate, given by Equation (19), and their CDF are shown in Figure 
5. In general, there is no obvious evidence for a violation of the regression assumptions. However, there 
appears to be a possible trend in the residuals with FWHM. To test this we calculated the linear correlation 
coefficient between the residuals and \ogFWHM . This correlation was significant at only « 1.1 a. Therefore 
there is no significant evidence for a correlation in the residuals with FWHM, justifying the assumption 
Mbh oc L^'^^gFW H M'^ . Rank correlation tests gave similar results. 

Very similar results were found by Vestergaard & Peterson (2006). However, our work differs from their 's 
in that Vestergaard & Peterson (2006) assumed that the C IV and H/3 BLRs have the same dependence on 
Ljjv, R oc , and only fit the constant term. In contrast, we fit both the constant term and the coefficient 
for the dependence of line width on luminosity. 

One can use Equation (9) to combine the mass estimates based on the Mbh-L with the C IV broad 
line mass estimate. Combining the two relationships, and taking a{z) oo, we find 

/i = -30.53 + 0.54/ + 1.40 log i^VFifM. (20) 

The intrinsic uncertainty in /x is reduced to ct^ « 0.33 dex, an improvement of « 18% over the broad line 
mass estimates. However, as mentioned in § 6.1, there is considerable systematic uncertainty on the behavior 
of the Mbh-L relationship outside of the range probed by the reverberation mapping sample, and thus it 
may be safer to just use rhciv- 

There has been some concern over the effectiveness of using the C IV emission line for estimating quasar 
black hole masses (Shemmer et al. 2004; Baskin & Laor 2005). If one ignores the C IV FWHM, and only 
uses the continuum luminosity to estimate Mbh, then the intrinsic scatter is cri « 0.6 dex (cf. § 6.1). 
However, the intrinsic uncertainty in an estimate of Mbh using the C IV R-L relationship and FWHM is 
o'civ ~ 0.4 dex, a improvement over the Mbh-L relationship of 1/3. This improvement is significant, 
and therefore we conclude that the C IV line may be used for estimating SMBH masses. 

We can estimate by comparing the C IV FWHM for the RMS spectrum with that obtained from 
a single-epoch spectrum. To do this we use the sources from Peterson et al. (2004) with C IV FWHM 
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measured from the RMS spectrum, FWHMjims, excluding Fairall 9. We omitted Fairall 9 because its 
FWHMrms was considered to be unreliable by Peterson et al. (2004). The remaining five data points are 
consistent with a 1:1 relationship between FWHMrms and FWHMse- Estimating the intrinsic scatter 
using these five data points is diflacult, particularly because the measurement errors are large. In fact, 
the observed scatter about FWHMjims ~ FWHMse is consistent with entirely being the result of the 
measurement errors. This suggests that is small, and it is likely that (t„ < 0.1 dex. 

If we assume that the only sources for scatter in the broad line estimates of m are from the R L 
relationship and from using the single-epoch line width to estimate V, then a value of ay « 0.1 dex implies 
that the intrinsic scatter about the R-L relationship for C IV is ar ~ 0.35 dex. This is about a factor of two 
larger than the value (7^ « 0.18 dex found for H/? by Peterson et al. (2004), where we have converted the 
value ar « 40% to dex. This larger scatter in the C IV R-L relationship implies that the H/3 mass estimates 
are more efficient than the C IV ones. However, Vestergaard & Peterson (2006) found that the H/3 mass 
estimates have a statistical scatter of cth/j ~ 0.43 dex, and thus are just as accurate as the C IV-based ones. 
This larger value of anp implies that cither there is significant uncertainty in using the SES H/3 line width 
to estimate the BLR velocity dispersion, a^^ ~ 0.36, or that there are additional sources of uncertainty in 
the broad line mass estimates beyond the scatter in the R-L relationship and the uncertainty in using the 
SES line width instead of the variable component line width. If there arc additional sources of scatter in the 
SES mass estimates, then ar ~ 0.35 dex represents an upper bound on the C IV R-L relationship scatter. 

Equation (7) implies that the scatter in the R-L relationship is the result of variations in the BLR 
ionization parameter, gas density, and average ionizing photon energy. While there are likely other sources 
of scatter in the R-L relationship, such as inclination, variability, and conversion between L\ and Lion, 
differences in BLR properties certainly contribute to this scatter and possibly dominate it. The fact that 
o'^^^ ^ 1.6o-^^ suggests that the magnitude of the dispersion in these properties is larger for the C IV 

1 /2 

emitting region. In addition, the regressions of (Peterson et al. 2004) are consistent with i?H/3 oc ^1450, 
similar to the result found here and by Vestergaard & Peterson (2006) for C IV. 



7. COMPARING THE H/3- AND C IV-BASED MASS ESTIMATES 

It is worth comparing our estimates of m based on single-epoch C IV to those based on combining 
single-epoch H/3 with the R^i^-L relationship. To do this, we took our sample of 100 sources that have data 
for both C IV and H/3 and calculated estimates of m from each. We compare the mass estimates obtained 
from both emission lines, and find that the two give mass estimates that are consistent so long as one uses 
the UV continuum to estimate the BLR size. 

For the H/? emission line, we estimate the BLR size from both the optical and UV Ru^-L relationships. 

The BLR size, Ru/S, is estimated as 

-Ro(-,;:,^^y. (21) 



10 It - dy 'V lO^^ergs s 

For Ro and a in the optical Rui3-L relationship, we use the luminosity at 5100A and the average of the BCES 
bisector and FITEXY (Press et al. 1992) fits of Kaspi et al. (2005), (i?o, a) = (2.23±0.21,0.69±0.05). For the 

o 1/2 

UV Rhp^L relationship, we use the luminosity at 1450A and {Rf),a) = (2.38,0.5). We assume i?H/3 oc ^^450 
to allow more direct comparison with the C IV-based mass estimates, and because the fits of Kaspi et 
al. (2005) are consistent with this form. Using the averaged values of R for the Balmer lines as listed by 
Kaspi et al. (2005), we recalculated the R-L1450 relationship for H/3 with a fixed at a = 0.5, and found 
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Ro = 2.38 ± 0.25. We then use Ru/} with H/3 FWHM to calculate toh/3 from Equation (1), after converting 
the H/3 FWHM to a velocity dispersion assuming FWHM/V = 2, based on the average FWHM/a* from 
Peterson et al. (2004). 

Based on the results of § 6.2, we calculate estimates of Mbh from the C IV line for these sources as 
"rhciv = —21.92 + O.5Z1549 + 2vciv- This form was found using the same method as described in § 6.2, but 
with 9i fixed at Oi = 1/2. In Table 5 we report the broad line mass estimates based on the Hp and C IV 
lines for these sources, as well as a weighted average of the two. 

In Figure 6 wc show mu/? vs mciv for both the optical and UV Rup-L relationship. We also show 
the 99% (2.6(7) confidence region for the parameters of a BCES bisector fit to riiciv = A + Brhu/B- As 
can be seen, the values of toh/3 calculated using the UV Rh0-L relationship are more consistent with the C 
IV-based ones. In addition, a pure 1:1 relationship falls outside of the 99% confidence region on the BCES 
bisector parameters for the optical-based mug values, but it is contained within the 99% confidence region 
for the UV-based estimates. In fact, the BCES bisector fit for the UV-based mass estimates differs from a 
1:1 relationship at a significance level of only 1.90-. Our results are in agreement with Warner et al. (2003) 
and Dietrich & Hamann (2004), who also compared the C IV- and H/3-based mass estimates and found them 
to be consistent. 

That the C IV- and II/3-based mass estimates are consistent is in contrast to the conclusions of Baskin 
& Laor (2005) and Shemmer et al. (2004). Baskin & Laor (2005) suggested that C IV-based mass estimates 
may be problematic because they find that \og{FWHMciv/FWHMii0) is significantly anti-correlated with 
log FW H Mufj . However this is expected even when vup oc ?;civ , as an anti-correlation between the ra- 
tio of the line widths and the H/3 FWHM is expected by construction; i.e., an anti-correlation between 
log FWHMcw — log FWHMuf} and log FWHMu^ is expected so long as the deviations of log FWHMcw 
and log FWHMu/3 from there respective means are not strongly correlated. In addition, none of the sources 
in the Baskin & Laor (2005) sample have H/3 broader than C IV for FWHM^^ < 4000 km s~^. To a lesser 
extant, a similar efi'ect is seen in Figure 7 of Shemmer et al. (2004) and Figure 4 of Warner et al. (2003). 
In Figure 6 we compare the C IV and H/3 FWHM for our sources. We find that most of the sources with 
FWHMYifj < 2000 km s^^ tend to have broader C IV lines. While the line widths certainly do not have a 
1:1 relationship, the BCES bisector fit found that on average FW H MB_fj is approximately proportional to 
FWHMciv, FWHMciv oc FWHM^P^° °^. Baskin & Laor (2005) and Shemmer et al. (2004) did not 
perform a regression analysis, so we are unable to make a more quantitative comparison. 

If the line widths are set by the virial relationship, we would expect FWHMciv oc FWHMuis. While 
the divergence from FWHMciv oc FW H M-ap is likely real (3.5<j significance), it is small, especially when 
compared with the intrinsic scatter in the FWHM plot. Furthermore, this divergence from FWHMciv oc 
FWHMufi does not appear to have a significant effect on the C IV mass estimates, as the masses inferred 
from the two emission lines differ only at the level of 1.9cr. As can by seen from the plot in Figure 6, any 
systematic difference between the C IV- and H/3-based mass estimates is small compared to the statistical 
scatter in m. 

Vcstcrgaard & Peterson (2006) performed a rcanalysis of the Baskin & Laor (2005) sample and concluded 
that the poor correlation between the C IV FWHM and the H/3 FWHM seen by Baskin & Laor (2005) 
was due to the inclusion of a large number of Narrow Line Seyfert Is (NLSl) and the subtraction of a C IV 
narrow component. Vestergaard (2004) noted that for these sources the C IV line profile is unlikely to be 
suitable for estimating Mbh, as there may be a strong component from an outflowing wind. Upon removing 
the NLSls and the sources with poor lUE data, Vestergaard & Peterson (2006) found the correlation between 
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C IV FWHM and H/3 FWHM to be significantly better. 

In addition, Baskin & Laor (2005) compared the correlation of C IV EW with L/Ledd oc L/Mbh- A 
similar investigation was performed by Shemmer et al. (2004), where they compared the mctallicity indicator 
N V/C IV with L/Ledd- Both authors found a stronger correlation when Mbh was estimated from the H/3 
line as compared to C IV, and suggested that the C IV line may give a less efficient and possibly biased 
estimate of Mbh. However, as noted in § 2.3, correlation coefficients inferred from the estimated black 
hole masses must be interpreted with caution. In particular, the SES estimate of m, m, merely defines the 
centroid of the probability density of m|/,v,z, and is in general not the actual m for an object. Therefore, 
this difference in the correlation coefficient between the H/3- and C IV-based mass estimates may be just the 
result of random sampling, and it is unclear whether one can say that the correlation coefficients inferred 
from the two different mass estimates are inconsistent with those mass estimates being drawn from the same 
parent distribution. 

The most direct test of the effectiveness of using C IV to estimate Mbh is found by comparing those 
values of rh estimated using C IV with the actual reverberation mapping values of m. This was done in § 6.2, 
where inclusion of the C IV FWHM resulted in a reduction of the error in m of ^ 1/3. In addition. Figure 
5 compares the C IV-based estimate of m with the reverberation mapping values. As can be seen, assuming 
Mbh oc FWHM^jy is consistent with the reverberation mapping sample, and therefore there is no reason 
to assume that the C IV line gives a biased estimate of Mbh- If inciv does give a biased estimate of m, 
then this bias is likely negligible compared to the variance in rhciv- 

8. SUMMARY 

In this work we have undertaken a statistical investigation of the best method to estimate SMBH mass 
based on the single-epoch C IV line and AGN continuum at 1549 A. The main conclusions are : 

1. Estimating AGN black hole masses from a single-epoch spectrum is a considerably different problem 
than in the reverberation-mapping case. Because of this, one is not estimating Mbh directly, but 
calculating a probability distribution of Mbh given the observed luminosity, line width, and redshift, 
p{m\l, V, z) oc p{v\l, m, z)p{l\m, z)p{m\z). 

2. Combining the information in L from both an intrinsic Mbh-L correlation and the R-L relationship 
results in improved black hole mass estimates. However, bccaiise of the current systematic uncertainties 
in the Mbh-L relationship, estimates based on its inclusion should be viewed with caution. In addi- 
tion, incorporating information from the intrinsic distribution of Mbh also results in more accurate 
estimates, on average. 

3. The distribution of Mbh inferred from the broad line mass estimates, msL, is broader than the intrinsic 

distribution of Mbh, as to^l are estimates of m contaminated by 'measurement error.' In addition, 
it is necessary to estimate or assume p(l\m, z) in order to estimate a survey's completeness in Mbh- 

4. The best estimates of Mbh based on the C IV emission line are obtained using the UV contimmm 
luminosity and FWHM. We find evidence that the C IV line shift, centroid, EW, and spectral slope 
of the UV continuum do not contribute any additional information to estimating Mbh for a given 
Li549 and FWHMciv- 
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5. Using the reverberation mapping sample, we find an Mbh-L relationship of the form 

/ = 35.72(±1.67) + 1.17(±0.22)m. (22) 

The intrinsic scatter about this relationship is cr; w 0.61 dex. Combining mass estimates based on this 
relationship with the broad line mass estimates results in a reduction in the statistical error of ~ 18%. 

6. We estimate a C IV R-L relationship of the form R a L^^^^^'^"^ . This is consistent with an R-L 

relationship of the form R cx L}/"^ predicted from simple photoionization physics or if the BLR size is 
set by the dust sublimation radius. The scatter about the C IV R-L relationship inferred from the 
reverberation mapping estimates of Mbh is Ur ~ 0.35 dex. 

7. A broad line estimate of m based on the C IV emission line may be calculated as 

logMBH/MQ = rhciv = -17.82(±2.99) + 0.41(±0.07) logAL;^(1549l) + 2logFWHMciv (23) 

Here, XLx is in imits of ergs and FWHM is in units of km s^^. The correlation between regression 
coefficients is Corr{mQ^^ ,9i) = —0.9996. The intrinsic scatter in m about Equation (23) is « 0.40 
dex. 

1 /2 

8. The C IV- and H/3-based mass estimates are consistent if one assumes R oc L^y for both emission 
lines. The two emission lines give estimates of Mbh with comparable accuracy. 

9. We find evidence that the R-L relationships for C IV and II/3 are similar in their dependence on Ljjv- 

1 /2 

For both, R oc LjjY is consistent with the available data. Also, the scatter in the C IV R-L relationship 

is a factor < 1.6 larger. This suggests that the C IV broad line region gas has a larger dispersion in 
its properties (e.g., density, ionization parameter, inclination) than the H/3 BLR gas. 
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Fig. 1. Results of the simulations described in § 2.2, illustrating the cfFectivcness of the mass estimate, 
/X, that combines the broad line mass estimates with the distribution of luminosities at a given black hole 
mass and the intrinsic distribution of Mbh- The upper panel shows the distribution of the residuals when 
using fi (dashed line) and the broad line mass estimate, rfiBL (solid line). The bottom panel compares the 
intrinsic distribution of m (thick solid line) with that inferred from the distribution of /x (dashed line) and 
niBL (thin solid line). 
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Fig. 2a. — New optical spectra, shown in the observed frame. The fluxes are in units of 10~^^ ergs cm~^ 
sec"'^ A~^. The absorption features at ^ 6875A and ~ 7600A are the A- and B-band atmospheric absorption 
hues. Note the strong iron emission in Q1340-0038. 
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Fig. 2b.— 
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Fig. 3. — Plots showing m as a function of the measured quantities for the sample with reverberation 
mapping data. Three outliers in the m vs a plot have been removed so as to make the structure in the plot 
easier to see. 
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Fig. 4. — The residuals I ~ I as a function of m for the Mbh^L relationship and the empirical CDF of 
the standardized residuals, {I — 1)/(Ti. The dashed lines define the 95% pointwise confidence interval of the 
empirical CDF, and the smooth solid line is the CDF of the standard normal density. 
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Fig. 5. — The residuals m — rhciv for the C IV broad Une mass estimate, shown as a function of mciv, I, 
and FWHM. Also shown is the empirical CDF of the standardized residuals, (to — rhciv)/o'civ ■ The 
dashed lines define the 95% pointwise confidence interval of the empirical CDF, and the smooth solid line is 
the CDF of the standard normal density. 
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Fig. 6. Plots comparing the H/3-based mass estimates with the C IV-bascd ones. The values of toh/3 for 
the top left plot are calculated from the optical R-L relationship, while the values of mn/S for the top right 

1/2 

plot assume i?H/3 L-^'^^q. The error bars denote the propagated measurement errors in rh resulting from the 
measurement errors in the line widths and continuum luminosities. The solid lines are a 1:1 relationship and 
the dashed lines are the BCES bisector fits. The bottom left plot shows the 99% (2.6cr) confidence regions 
for the BCES bisector fits comparing rhciv with rriH/S- The solid contour correspond to rhuis calculated 
from the optical R}ip-L relationship, and the dashed contour when toh/3 is calculated from the UV one. 
The crosses mark the best-fit values, and the triangle marks the values expected for a 1:1 relationship. The 
bottom right plot shows a comparison between the H/3 and C IV FWHM. The solid line is the best fit 
FWHMu/s oc FWHMciv relationship, and the dashed line is the BCES bisector fit. 
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Table 1. Object List 



Object 


RA (J2000) 


DEC (J2000) 


Redshift 


log ( Mb H 




Rad. Type*" 


Instrument 


MRK 335 


00 


06 


19.5 


+20 12 10.5 


0.025 


6.41 


± 


0, 


.11 


Quiet 


FOS 


PG 0026+129 


00 


29 


13.7 


+13 16 03.8 


0.142 


7.85 


± 


0, 


.10 


Quiet 


FOS 


PG 0052+251 


00 


54 


52.1 


+25 25 39.3 


0.155 


7.82 


± 


0, 


.08 


Quiet 


FOS 


FAIRALL9 


01 


23 


45.7 


-58 48 21.8 


0.046 


7.66 


± 


0, 


.09 


Quiet 


FOS 


MRK 590 


02 


14 


33.6 


-00 46 00.1 


0.027 


6.93 


± 


0, 


.06 


Quiet 


lUE 


3C 120 


04 


33 


11.1 


+05 21 15.6 


0.033 


7.00 


± 


0, 


.21= 


Loud 


lUE 


ARK 120 


05 


16 


11.4 


-00 08 59.4 


0.033 


7.43 


± 


0, 


.05 


Quiet 


FOS 


MRK 79 


07 


42 


32.8 


+49 48 34.7 


0.0221 


6.97 


± 


0, 


,11 


Quiet 


lUE 


PG 0804+761 


08 


10 


58.6 


+76 02 42.0 


0.1000 


8.10 


± 


0, 


,05 


Quiet 


lUE 


MRK 110 


09 


25 


12.9 


+ 52 17 10.5 


0.0352 


6.65 


± 


0, 


,10 


Quiet 


lUE 


PG 0953+414 


09 


56 


52.4 


+41 15 23.0 


0.2341 


7.69 


± 


0, 


,09 


Quiet 


FOS 


NGC 3516 


11 


06 


47.5 


+72 34 06.9 


0.0088 


6.88 


± 


0, 


,14 


Quiet 


FOS 


NGC 3783 


11 


39 


01.7 


-37 44 18.9 


0.009 


6.73 


± 


0, 


,07 


Quiet 


FOS 


3C 273.0 


12 


29 


06.7 


+02 03 09.0 


0.1583 


8.20 


± 


0, 


,09 


Loud 


FOS 


PG 1307+085 


13 


09 


47.0 


+08 19 49.8 


0.1550 


7.90 


± 


0, 


12 


Quiet 


FOS 


MRK 279 


13 


53 


03.4 


+69 18 29.6 


0.0304 


6.80 


± 


0, 


26 


Quiet 


STIS 


NGC 5548 


14 


17 


59.5 


+ 25 08 12.4 


0.0171 


7.08 


± 


0, 


,01 


Quiet 


FOS 


PG 1426+015 


14 


29 


06.6 


+01 17 06.5 


0.0864 


8.37 


± 


0, 


,12 


Quiet 


lUE 


MRK 817 


14 


36 


22.1 


+ 58 47 39.4 


0.0314 


6.95 


± 


0, 


,06 


Quiet 


lUE 


PG 1613+658 


16 


13 


57.2 


+65 43 09.6 


0.1290 


7.70 


± 


0, 


20'^ 


Quiet 


lUE 


PG 1617+175 


16 


20 


11.3 


+ 17 24 27.7 


0.1124 


8.03 


± 


0, 


10 


Quiet 


lUE 


3C 390.3 


18 


42 


08.8 


+79 46 17.0 


0.0560 


7.71 


± 


0, 


,09 


Loud 


FOS 


MRK 509 


20 


44 


09.8 


-10 43 24.5 


0.0350 


7.41 


± 


0, 


,03 


Quiet 


FOS 


PG 2130+099 


21 


32 


27.8 


+ 10 08 19.5 


0.0629 


7.91 


± 


0, 


,05 


Quiet 


GHRS 


NGC 7469 


23 


03 


15.6 


+08 52 26.4 


0.0163 


6.34 


± 


0, 


,04 


Quiet 


FOS 



^Black hole masses are from Peterson et al. (2004), assuming / = 1. 
^The radio type is either radio loud or radio quiet. 

'^The standard errors in Mbh for 3C 120 and PG 1613+658 were estimated by averaging the upper and lower 
uncertainties listed by Peterson et al. (2004). 
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Table 2. Log of New Observations 



Object 


RA (J2000) 


DEC (J2000) 


Redshift 


Date 


Instrument 


Exp. Time (s) 


Q 0003+0146 


00 


05 


47.6 


+02 03 02.2 


0.234 


Oct 06, 2002 


Bok B&C 


1800 


Q 0017+0209 


00 


20 


25.1 


+02 26 25.3 


0.401 


Oct 06, 2002 


Bok B&C 


2700 


PKS 0044+030 


00 


47 


05.9 


+03 19 54.9 


0.624 


Oct 06, 2002 


Bok B&C 


1800 


Q 0100+0205 


01 


03 


13.0 


+02 21 10.4 


0.394 


Oct 06, 2002 


Bok B&C 


2700 


Q 0115+027 


01 


18 


18.5 


+02 58 05.9 


0.672 


Oct 06, 2002 


Bok B&C 


3600 


3C 57 


02 


01 


57.2 


-11 32 33.7 


0.669 


Oct 06, 2002 


Bok B&C 


1800 


PKS 0214+10 


02 


17 


07.7 


+ 11 04 09.6 


0.408 


Oct 06, 2002 


Bok B&C 


2700 


IR 0450-2958 


04 


52 


30.0 


-29 53 35.0 


0.286 


Oct 06, 2002 


Bok B&C 


207 


3C 232 


09 


58 


21.0 


+32 24 02.2 


0.533 


May 15, 2002 


Bok B&C 


1800 


4C 41.21 


10 


10 


27.5 


+41 32 39.1 


0.611 


May 15, 2002 


Bok B&C 


1800 


B2 1028+313 


10 


30 


59.1 


+31 02 56.0 


0.178 


May 15, 2002 


Bok B&C 


1800 


MC 1104+167 


11 


07 


15.0 


+ 16 28 02.4 


0.632 


May 15, 2002 


Bok B&C 


1800 


Q 1132-0302 


11 


35 


04.9 


-03 18 52.5 


0.237 


May 14, 2002 


Bok B&C 


1800 


Q 1150+497 


11 


53 


24.5 


+49 31 08.6 


0.334 


May 14, 2002 


Bok B&C 


1200 


TEX 1156+213 


11 


59 


26.2 


+21 06 56.2 


0.349 


May 15, 2002 


Bok B&C 


2100 


MRK 205 


12 


21 


44.0 


+75 18 38.1 


0.070 


May 15, 2002 


Bok B&C 


1800 


Q 1230+0947 


12 


33 


25.8 


+09 31 23.0 


0.420 


Mar 20, 2004 


Bok B&C 


1800 


Q 1317-0142 


13 


19 


50.3 


-01 58 04.6 


0.225 


May 14, 2002 


Bok B&C 


900 


PG 1333+176 


13 


36 


02.0 


+ 17 25 13.0 


0.554 


May 15, 2002 


Bok B&C 


2100 


Q 1340-0038 


13 


42 


51.6 


-00 53 46.0 


0.326 


May 14, 2002 


Bok B&C 


1200 


Q 1401+0951 


14 


04 


10.6 


+09 37 45.5 


0.441 


May 15, 2002 


Bok B&C 


1800 


Q 1425+2003 


14 


27 


25.0 


+ 19 49 52.3 


0.111 


May 14, 2002 


Bok B&C 


900 


PKS 1451-375 


14 


54 


27.4 


-37 47 34.2 


0.314 


Jul 17, 2004 


IMACS 


2700 


PG 1538+477 


15 


39 


34.8 


+47 35 31.6 


0.770 


May 14, 2002 


Bok B&C 


1500 


PKS 2128-12 


21 


31 


35.4 


-12 07 05.5 


0.501 


Oct. 06, 2002 


Bok B&C 


1800 


PKS 2135-147 


21 


37 


45.2 


-14 32 55.8 


0.200 


May 15, 2002 


Bok B&C 


900 


Q 2141+175 


21 


43 


35.6 


+17 43 49.1 


0.213 


May 15, 2002 


Bok B&C 


1200 


PKS 2243-123 


22 


46 


18.2 


-12 06 51.2 


0.630 


Oct 06, 2002 


Bok B&C 


1800 


PKS 2352-342 


23 


55 


25.6 


-33 57 55.8 


0.706 


Jul 17, 2004 


IMACS 


1800 
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Table 3. Continuum and Fe Emission Fitting Windows 





1 




2 




3 




4 




5 




6 


uv 


AA1350- 


-1365 


AA1427- 


-1500 


AA1760- 


-1860 


AA1950- 


-2300 


AA2470- 


-2755 


AA2855-3010 


Optical 


AA3535- 


-3700 


AA4100- 


-4200 


AA4400- 


-4700 


AA5100- 


-6200 


AA6800- 


-7534 
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Tabk' 4. Coutinuum and IV Emission Lin(' ParanK^tcrs 



Object 


log ALa(1549 A) 


FWHM 




EW 


Ad 




a 




ergs cm^^ scc^"^ 


1000 km s"^ 


1000 km 


A 


1000 km 


1000 km 




MRK 335 


44.38 


2.094 


1.856 


71.82 


0.368 


0.141 


-0.45 


PG0026+12 


45.19 


1.971 


5.263 


52.45 


0.054 


-0.99 


-0.85 


0052+2509 


45.39 


5.762 


4.544 


120.4 


-0.63 


-0.51 


-0.29 


FAIRALL9 


44.39 


2.686 


2.633 


193.0 


-0.19 


-0.02 


-0.72 


MRK 590 


43.35 


2.846 


1.818 


185.2 


-0.72 


-1.21 


-0.80 


3C 120 


44.56 


5.040 


2.813 


148.3 


-1.34 


-0.87 


-0.77 


ARK 120 


43.91 


4.127 


3.124 


70.0 


-0.00 


0.647 


-0.95 


MRK 79 


43.68 


4.615 


2.891 


112.4 


-0.41 


-0.68 


-0.73 


PG 0804+761 


45.26 


3.608 


3.307 


62.63 


-0.41 


-1.20 


-0.79 


MRK 110 


43.77 


3.468 


2.117 


185.0 


0.819 


0.437 


-0.79 


PG0953+414 


45.49 


3.114 


4.339 


95.38 


-1.20 


0.210 


-0.25 


NGC 3516 


42.31 


3.512 


2.480 


327.9 


0.054 


-0.47 


-2.05 


NGC 3783 


43.23 


2.376 


2.527 


110.8 


0.054 


-0.21 


-0.51 


3C 273.0 


46.28 


3.659 


3.285 


31.24 


-1.70 


-0.64 


-0.57 


1307+0835 


45.08 


4.054 


3.395 


96.97 


0.431 


0.339 


-0.35 


MRK 279 


44.32 


3.057 


3.201 


110.2 


-0.85 


-1.00 


-0.29 


NGC 5548 


43.49 


3.490 


3.278 


202.2 


0.431 


0.214 


-0.69 


PG 1426+015 


45.21 


5.398 


4.422 


43.81 


-0.10 


-0.00 


-0.61 


MRK 817 


44.13 


4.192 


3.719 


66.52 


-0.10 


0.153 


-0.94 


PG 1613+658 


45.14 


5.869 


5.099 


81.43 


-0.10 


1.562 


-0.79 


PG 1617+175 


44.69 


7.842 


2.906 


50.60 


0.201 


0.947 


-0.75 


3C390.3 


43.81 


6.325 


5.113 


138.2 


-0.13 


-0.78 


1.62 


MRK509 


44.68 


4.313 


3.722 


128.3 


-0.95 


-0.13 


-0.77 


PG 2130+099 


44.60 


2.200 


1.785 


58.36 


-0.78 


-0.94 


-0.78 


NGC 7469 


43.29 


2.822 


2.420 


161.3 


-0.07 


0.028 


-1.75 
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Table 5. Black Hole Mass Estimates 



RA (J2000) 


DEC (J2000) 


z 


logM^^/Mo" 




log M^^J/Mq 


00 


05 


47.6 








7 RI _|_ n AC\ 


7 /17 -1- n OQ 


00 


05 


59.2 


11 ft no AQ 1 


450 




Q 1 Q _|_ n AC\ 
0. J.O ZLL U.IU 


8 97 -1- n OQ 
o.z 1 zl u. zy 


00 


06 


19.5 


_i_on 10 1 n 


n noc; 
u.u.^o 


O. (.^ i U.4o 


R on -\- (\ Ar\ 
o.yu i U.4U 


R 80 -1- n OQ 
o.oz i u.zy 


00 


20 


25.1 




401 


7 oo _|_ n 40 


8 'X'X -\- n An 

0.00 ZLL U.4U 


7 89 -1- n 9Q 

/ .oz zll u . zy 


00 


29 


13.7 


H) U,>.o 


142 


7 -1- n /I 


7 1 7 -1- n /I n 

/ . 1 / i U.4U 


7 9^ -1- n 9Q 
( .ZO i u.zy 


00 


47 


05.9 




624 


Q n7 -1- n A'\ 


8 fi7 -1- n An 

0.0/ ^ U.4U 


a afi -1- n 9q 
o.oD zll u.zy 


00 


52 


02.4 


^ni m OQ 


2 270 


Q A'7 -\- r\ AA 


Q /1 1 _i_ n /I n 

y,41 i U.4U 


Q /I /I -1- n 9Q 

y.44 i u.zy 


00 


52 


33.7 


u_ni zin An fi 


2 307 


Q 91^ -1- n Af, 
y . ZO ^ u.4:D 


n -|- n An 
y.oo ^ U.4U 


Q A9 -1- n "^n 
y . 'rt z ^ u . u 


00 


54 


52.1 


-f^zo zo oy.o 


155 


7 CO _i_ n A'i 
1 . oy ^ u.4:o 


8 '^A -1- n An 

0.04 ZLL U.4U 


a 1 -U n 9Q 
. -LO zll u . zy 


01 


03 


13.0 


^n9 91 inA 


394 


s 47 -1- n A'i 


Q 1 c _i_ n An 

. -LD ZLL U.4U 


a '^n -1- n 9q 
.OU zll u . zy 


01 


18 


18.5 


^n9 Q 


672 


Q OQ _|_ ri A A 
O.Oo ^ U.44 


7 CI _i_ n An 

/ . D 1 ZlL U.4U 


7 Qf, -U n 9Q 

/ . yo zll u . zy 


01 


23 


45.7 


S ,1 S 91 S 
-Oo Zl.o 


047 


T AC\ 4- n A'i 
1 .4U i U.40 


7 1 c -1- n /in 

/ , ID i U.4U 


7 97 -1- n 9Q 
/ .z / i u.zy 


01 


26 


42.8 


^9^ ni 


2 370 


Q 9zl -1- n 


n 9c, _|_ n An 


Q 9'^ -1- n "^n 
y . ZO zll u . OU 


01 


57 


35.0 


-\- i ^ 4iZ 4iO.Z 


Z.OOO 


Q 1 n -1- n /I /I 
y . i u i u . 44 


Q 9n -1- n /1 1 
y .zu i U.41 


Q 1 c _l_ n "^n 
y . 10 i u.ou 


02 


01 


57.2 


-11 32 33 7 


669 


c ni -1- n A7 

. U-L ZLL U.4 ( 


Q Qc _|_ n An 

O.oO ZLL U.4U 


a AQ -1- n "^n 
o.4y ^ u.ou 


02 


17 


07.7 


^11 ri/i riQ fi 


408 


Q CO _|_ n /I ^ 
.oy i u.40 


a 9n -1- n /I n 

o.ZU i U.4U 


a /1 9 -1- n 90 
0.4Z i u.zy 


02 


59 


05.6 


^nn 11 9 1 Q 


3 366 


Q or; _|_ fi A'^ 

. y zll u .40 


Q n'^ -1- n An 

y.uo ^ u.4:U 


a QQ -U n 9Q 
. y y ^ u.zy 


03 


04 


49.9 


nn ns i /i 




Q on -1- n /I 
o.yu i u.40 


Q en _|_ n /I n 
o.oy i U.4U 


a 7Q -1- n 9Q 
0. ( y i u.zy 


03 


51 


28.6 


1 A 9Q flQ 1 
zy uy.j- 


616 


8 Q9 -1- n A"^ 
o.yz zll U.40 


Q n9 -1- n An 
y.uz ^ U.4U 


a Qa -1- n 9q 
. yo ^ u . zy 


04 


05 


34.0 


1 ns 1/11 

-lO Uo 14.1 


n 1^71 
y.1.0 1 1 


.OD i u.40 


a ni^ -1- n /I n 

O.UO i U.4U 


a 9Q -1- n 9Q 
o.zy i u.zy 


04 


07 


48.4 


19 11 '^fi n 

-IZ 11 ou.u 


574 


8 91^ -1- n A"^ 

. ZO ^ U.40 


8 Ai -1- n An 

0.4± ^ U.4U 


a '^A -1- n 9Q 
. 04 zll u . zy 


04 


17 


16.8 


nf^ /I Q 
-uo oo 4io.y 


n Ts 1 


Q C\A 4- C\ 
y.U4 i U.40 


Q QR _i_ n /I n 
o.yo i U.4U 


Q nn -1- n "^n 
y.uu i u.ou 


04 


41 


17.3 


A'i 1 7 
lO 40./ 


593 


8 ns -1- n A*? 

O.Uo ZLL U.40 


7 8A -1- n An 

/ .o4 ^ U.4U 


7 Qf^ -1- n 9Q 
1 . yo ^ u . zy 


04 


52 


30.0 


9Q "i^ n 
-zy oo oo.u 


U.ZoD 


7 ^7 -1- n /I "i 

/ .0 1 i u.40 


7 8fi -1- n /I n 

( .oD i U.4U 


7 fi/i -1- n 9Q 
/ .D4 i U.zy 


04 


56 


08.9 


91 f^Q nQ A 

-zi oy uy.4 


534 


R 7*^ -1- n A^ 

. i ^ U.4:0 


8 AQ -1- n An 

0.4y ZLL U.4:U 


8 fin -1- n 9Q 
.Du ^ u . zy 


07 


45 


41.7 


-i-^1 AO fif^ 7 
-f-Ol tZ oo. 1 


462 


8 97 -1- n A*^ 

0. Z / ^ U.4:0 


8 ^A -1- n An 

0.04: ZLL U.4:U 


8 A1 -I- n 9Q 
.4 J. zii U.zy 


08 


40 


47.6 


11 Q 1 9*^ 7 

-f-lO IZ ZO. 1 


684 


8 nfi -1- n d.*^ 

O.UD ZlL U.^O 


Q nq _|_ fi ^fi 


8 07 -1- 9Q 
o.u ( zc u.zy 


08 


53 


34.2 


_i_A^ AQ ni 


513 


7 84 -1- n d.'^ 




8 1 c; -1- OQ 
0. J-O zc u. zy 


09 


19 


57.7 


_i_c^i nfi inn 
-f-oi uo lU.U 


u.ooo 




8 c;8 -1- n /in 


a /i'^ -1- n OQ 
0.40 It U.zy 


09 


50 


48.4 


ion Oft K-t n 
-f-oy ZO oi.u 


206 


0.00 ZlL u.^o 


7 88 -1- n 4n 
1 .00 ^ u.^u 


8 1 8 -1- OQ 
0. io zc U.zy 


09 


56 


52.4 


_i_A lit; OQ n 

-f-^l 10 ZO.U 


n OQO 
u.^oy 


7 7n -\- (\ A"^ 


7 8^ -1- n /in 


7 78 -1- n OQ 
( . ( i U.zy 


09 


58 


21.0 


loo OA no 

T^OZ Z'i UZ.Z 


533 




8 87 -1- n 40 
0.0 ( ZlL u.^u 


8 90 -1- "^0 
O.ZZ Zll U.OU 


10 


04 


02.6 


_1_0S2 KK "iK n 
-f-ZO oo «50.U 


n Qoo 
u.o.^y 


R 8R -1- n /l'^ 


8 ot^ -1- n /in 

o.ZO i U.4U 


7 R1 -1- n OQ 
1 .Dl i U.zy 


10 


04 


20.1 


inc; 1 o nn n 

-f-UO lo uu.u 


161 




7 40 -1- 40 


7 4Q -1- OQ 
1 .^y zii U.zy 


10 


10 


27.5 


_i_A1 ^0 39 1 

T^'il oz oy. 1 


0.611 


7 qq _l_ n 4*^ 


0.00 Lf.'iO 


a 90 -1- n 9Q 


10 


30 


59.1 


+31 02 56.0 


0.178 


7.97 ± 0.44 


7.59 ± 0.40 


7.76 ±0.30 


10 


51 


51.5 


+00 51 18.1 


0.357 


7.88 ± 0.43 


8.18 ± 0.40 


8.04 ±0.29 


11 


04 


13.9 


+76 58 58.2 


0.311 


8.44 ± 0.43 


7.92 ± 0.40 


8.16 ±0.29 


11 


06 


31.8 


+00 52 53.4 


0.425 


8.67 ±0.43 


8.43 ± 0.40 


8.54 ±0.29 


11 


06 


33.5 


-18 21 24.0 


2.319 


8.89 ± 0.45 


9.65 ± 0.40 


9.31 ± 0.30 


11 


07 


15.0 


+16 28 02.4 


0.632 


8.04 ±0.45 


8.65 ± 0.40 


8.38 ±0.30 


11 


18 


30.3 


+40 25 55.0 


0.154 


7.82 ± 0.43 


7.93 ± 0.40 


7.88 ±0.29 


11 


19 


08.7 


+21 19 18.0 


0.176 


7.55 ± 0.43 


8.13 ± 0.40 


7.87 ±0.29 


11 


24 


39.2 


+42 01 45.2 


0.234 


7.61 ± 0.43 


7.47 ± 0.40 


7.53 ± 0.29 


11 


35 


04.9 


-03 18 52.5 


0.237 


7.07 ± 0.43 


7.87 ± 0.40 


7.50 ± 0.29 


11 


39 


57.1 


+65 47 49.4 


0.652 


7.60 ± 0.43 


8.22 ± 0.40 


7.93 ± 0.29 


11 


41 


21.7 


+01 48 03.3 


0.383 


7.20 ± 0.45 


8.23 ± 0.40 


7.77 ± 0.30 


11 


47 


18.0 


-01 32 07.7 


0.382 


8.08 ± 0.43 


6.68 ± 0.43 


7.37 ± 0.30 


11 


53 


24.5 


+49 31 08.6 


0.334 


7.68 ± 0.43 


8.02 ± 0.40 


7.86 ± 0.29 


11 


58 


39.9 


+62 54 28.1 


0.594 


7.45 ± 0.51 


8.32 ± 0.40 


7.99 ± 0.32 


11 


59 


26.2 


+21 06 56.2 


0.349 


8.69 ± 0.43 


8.03 ± 0.40 


8.33 ± 0.29 


12 


04 


42.2 


+27 54 12.0 


0.165 


8.24 ± 0.43 


7.26 ± 0.40 


7.71 ± 0.29 


12 


14 


17.7 


+ 14 03 12.3 


0.080 


6.72 ± 0.43 


6.95 ± 0.40 


6.84 ± 0.29 


12 


19 


20.9 


+06 38 38.4 


0.334 


8.00 ± 0.43 


7.95 ± 0.40 


7.97 ± 0.29 
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Table 5 — Continued 



RA (J2000) 










log Mflf^/Mg 




12 


21 


44.0 


+ 75 18 38.1 


0, 


.070 


7.89 ± 


.43 


7.59 ± 


.40 


7.73 ± 0.29 


12 


31 


20.6 


+07 25 52.8 


2, 


.391 


9.14 ± 


.49 


9.19 ± 


.40 


9.17 ± 0.31 


12 


33 


25.8 


+09 31 23.0 


0, 


.420 


7.98 ± 


,43 


8.18 ± 


.40 


8.09 ± 0.29 


13 


01 


12.9 


+59 02 06.9 


0, 


.472 


7.91 ± 


,43 


9.01 ± 


.40 


8.50 ± 0.29 


13 


05 


33.0 


-10 33 20.4 


0, 


.286 


8.19 ± 


,43 


8.38 ± 


.40 


8.29 ± 0.29 


13 


09 


47.0 


+08 19 49.8 


0, 


.155 


8.54 ± 


,43 


7.87 ± 


.40 


8.18 ± 0.29 


13 


12 


17.7 


+35 15 21.0 


0, 


.184 


7.62 ± 


,43 


7.43 ± 


.40 


7.52 ± 0.29 


13 


19 


50.3 


-01 58 04.6 


0, 


.225 


7.05 ± 


,43 


7.23 ± 


.40 


7.15 ± 0.29 


13 


23 


49.5 


+65 41 48.0 


0, 


.168 


7.11 ± 


,43 


7.60 ± 


.40 


7.37 ± 0.29 


13 


36 


02.0 


+ 17 25 13.0 


0, 


.554 


7.72 ± 


,43 


8.79 ± 


.40 


8.29 ± 0.29 


13 


42 


51.6 


-00 53 46.0 


0, 


.326 


8.27 ± 


,43 


8.20 ± 


,40 


8.23 ± 0.29 


13 


57 


04.5 


+ 19 19 06.6 


0, 


.719 


7.12 ± 


,44 


7.99 ± 


,40 


7.60 ± 0.30 


14 


04 


10.6 


+09 37 45.5 


0, 


.441 


9.06 ± 


,45 


8.49 ± 


,41 


8.74 ± 0.30 


14 


05 


16.2 


+25 55 33.6 


0, 


,164 


7.35 ± 


,43 


8.19 ± 


,40 


7.80 ± 0.29 


14 


17 


00.9 


+44 56 06.0 


0, 


,114 


7.61 ± 


,43 


7.68 ± 


,40 


7.65 ± 0.29 


14 


19 


03.9 


-13 10 45.0 


0, 


,129 


7.56 ± 


,43 


7.76 ± 


.40 


7.67 ± 0.29 


14 


27 


25.0 


+19 49 52.3 


0, 


,111 


8.03 ± 0, 


.43 


8.30 ± 


.40 


8.17 ± 0.29 


14 


27 


35.7 


+26 32 15.0 


0, 


,366 


7.85 ± 0, 


.43 


8.80 ± 0, 


.40 


8.36 ± 0.29 


14 


29 


43.1 


+47 47 26.0 


0, 


,221 


7.65 ± 0, 


.43 


7.64 ± 0, 


.40 


7.64 ± 0.29 


14 


42 


07.5 


+35 26 22.9 


0, 


,077 


7.01 ± 0, 


.43 


7.26 ± 0, 


.40 


7.15 ±0.29 


14 


46 


45.9 


+40 35 07.1 


0, 


,267 


7.38 ± 0, 


.43 


8.42 ± 0, 


.40 


7.94 ± 0.29 


14 


54 


27.4 


-37 47 34.2 


0, 


,314 


7.62 ± 0, 


.43 


7.93 ± 


.40 


7.79 ±0.29 


15 


14 


43.5 


+36 50 51.0 


0, 


,371 


7.84 ± 0, 


.44 


8.28 ± 0, 


.40 


8.08 ±0.30 


15 


39 


34.8 


+47 35 31.6 


0, 


,770 


8.58 ± 0, 


.46 


9.01 ± 0, 


.40 


8.82 ± 0.30 


15 


45 


30.3 


+48 46 07.9 


0, 


,400 


7.73 ± 0, 


.43 


8.55 ± 0, 


.40 


8.17 ±0.29 


15 


47 


43.5 


+20 52 16.4 


0, 


,264 


8.39 ± 


.43 


8.21 ± 0, 


.40 


8.29 ±0.29 


16 


14 


13.2 


+26 04 16.2 


0, 


,131 


7.05 ± 0, 


.43 


7.78 ± 0, 


.40 


7.44 ± 0.29 


16 


20 


21.8 


+17 36 24.0 


0, 


,555 


8.22 ± 0, 


.43 


8.69 ± 0, 


.40 


8.47 ± 0.29 


16 


27 


56.1 


+55 22 31.0 


0, 


,133 


8.49 ± 0, 


.43 


7.89 ± 0, 


.40 


8.17 ± 0.29 


16 


42 


58.8 


+39 48 36.9 


0, 


,595 


8.01 ± 0, 


.43 


8.45 ± 0, 


.40 


8.24 ± 0.29 


17 


04 


41.3 


+60 44 30.0 


0, 


,371 


7.37 ±0, 


.44 


8.15 ± 0, 


.40 


7.80 ±0.30 


18 


21 


59.4 


+64 21 07.5 


0, 


,297 


8.85 ± 0, 


.43 


8.53 ± 0, 


.40 


8.68 ±0.29 


19 


27 


48.5 


+73 58 02.0 





,302 


8.13 ± 0, 


.43 


8.19 ± 0, 


.40 


8.16 ± 0.29 


20 


44 


09.8 


-10 43 24.5 





,035 


7.64 ±0, 


.43 


7.69 ± 0, 


.40 


7.66 ± 0.29 


21 


31 


35.4 


-12 07 05.5 


0, 


,501 


8.13 ± 


,43 


8.64 ± 


,40 


8.40 ± 0.29 


21 


37 


45.2 


-14 32 55.8 


0, 


,200 


8.86 ± 


,43 


7.99 ± 


,43 


8.42 ± 0.31 


21 


43 


35.6 


+ 17 43 49.1 


0, 


,213 


8.24 ± 


,43 


8.20 ± 


,41 


8.22 ± 0.30 


22 


03 


15.0 


+31 45 37.7 


0, 


,297 


8.08 ± 


,43 


8.76 ± 


,40 


8.45 ± 0.29 


22 


46 


18.2 


-12 06 51.2 


0, 


,630 


8.41 ± 


,43 


8.13 ± 


,40 


8.26 ± 0.29 


22 


54 


05.8 


-17 34 55.0 


0, 


,068 


8.15 ± 


,43 


7.56 ± 


,40 


7.83 ± 0.29 


22 


54 


10.4 


+ 11 36 38.9 


0, 


,323 


7.59 ± 


,43 


8.13 ± 


,40 


7.88 ± 0.29 


23 


03 


43.5 


-68 07 37.1 


0, 


,512 


7.62 ± 


,43 


7.77 ± 


,40 


7.70 ± 0.29 


23 


11 


17.8 


+ 10 08 16.2 


0, 


,432 


8.48 ± 


,43 


8.36 ± 


,40 


8.42 ± 0.29 


23 


46 


36.9 


+09 30 46.0 


0, 


,672 


7.94 ± 


,43 


7.99 ± 


,40 


7.97 ± 0.29 


23 


51 


56.0 


-01 09 13.7 


0, 


,174 


8.52 ± 0, 


.43 


8.09 ± 0, 


.40 


8.29 ± 0.29 


23 


55 


25.6 


-33 57 55.8 


0, 


,706 


6.51 ± 0, 


.44 


8.39 ± 0, 


.40 


7.53 ±0.30 



Note. — The Icr uncertainties include the contributions from measurement error in the Hne 
widths and continuum luminosities, and from the intrinsic uncertainty in 7%. 

^Thc H/3-bascd mass estimates arc for the UV Rufj—L relationship. The la errors were calcu- 
lated assuming cth/? — 0.43 dex (Vestergaard & Peterson 2006). 

^^BH ^ weighted average of the H/3- and C IV-based mass estimates. 



